Two-dimensional solitons with hidden and explicit vorticity in bimodal cubic-quintic 

media 
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We demonstrate that two-dimensional two-component bright solitons of an annular shape, carry- 
ing vorticities (m, ±m) in the components, may be stable in media with the cubic-quintic nonlin- 
earity, including the hidden-vorticity (HV) solitons of the type (m, — m), whose net vorticity is zero. 
Stability regions for the vortices of both (m, ±m) types are identified for m = 1, 2, and 3, by dint 
of the calculation of stability eigenvalues, and in direct simulations. A novel feature found in the 
study of the HV solitons is that their stability intervals never reach the (cutoff) point at which the 
bright vortex carries over into a dark one, hence dark HV solitons can never be stable, contrarily 
to the bright ones. In addition to the well-known symmetry-breaking {external) instability, which 
splits the ring soliton into a set of fragments flying away in tangential directions, we report two new 
scenarios of the development of weak instabilities specific to the HV solitons. One features charge 
flipping, with the two components exchanging the angular momentum and periodically reversing 
the sign of their spins. The composite soliton does not split in this case, therefore we identify such 
instability as an intrinsic one. Eventually, the soliton splits, as weak radiation loss drives it across 
the border of the ordinary strong (external) instability. Another scenario proceeds through separa- 
tion of the vortex cores in the two components, each individual core moving toward the outer edge 
of the annular soliton. After expulsion of the cores, there remains a zero-vorticity breather with 
persistent internal vibrations. 

PACS numbers: 42.65.Tg, 42.65.Jx, 42.65.Sf 



I. INTRODUCTION 

Optical vortices are two- or three-dimensional (2D or 
3D) dark or bright solitons with an embedded phase dis- 
location, which lends them a topological charge ( "spin" ) . 
As topologically nontrivial self-trapped states, the vor- 
tices have attracted a great deal of attention, which was 
additionally enhanced by a potential which the 2D ones 
may have as reconfigurable conduits for weak optical sig- 
nals P, • Vortex solitons of the bright type are of 
special interest in both respects. First, being relatively 
compact objects, they make it possible to realize sophis- 
ticated multi-vortex configurations. On the other hand, 
experimental creation of bright solitons is hampered by 
the fact that, in media with the simplest collapse- free 
nonlinearities (quadratic or saturable), they are subject 
to strong azimuthal instability, which splits them into a 
set of ordinary (zero-spin) solitons 0,O|. Nevertheless, 
as it was first observed in direct simulations of a model 
with the cubic-quintic (CQ) nonlinearity in Ref. 0, and 
later investigated in detail by more accurate methods 
(see a review in Ref. bright vortices with the spin 
to = 1 jSj, to — 2 ,3j, and m > 3 @ may be stable 
if the model features competing self-focusing and self- 
defocusing nonlinearities. Another example which cor- 
roborates this conclusion is a model combining quadratic 



and self-defocusing cubic nonlinearities [g. Later, stable 
two-component (vectorial) solitons, carrying the vortic- 
ity set (m, to), have been identified in a vectorial version 
of the CQ model @. 

It is relevant to mention that patterns of a similar type, 
in the form of "soliton necklaces" , i.e., ring-shaped chains 
of the fundamental (to = 0) solitons, were introduced 
and studied in the model with the Kerr (cubic) nonlin- 
earity The necklaces may be constructed with zero 
or non-integer angular momentum; however, they are not 
stationary objects, as they spread out and eventually dis- 
integrate. A possibility to (practically) stabilize necklace- 
like patterns is to introduce a vectorial interaction with 
a fundamental soliton, which helps to support vortex- 
mode, dipole-mode, and multipole vector solitons [Tl| . 
Additional stabilization of two-component necklaces in 
vectorial model allows to construct stationary necklace- 
ring vector solitons |l2j . having azimuthally modulated 
(necklace-type) components whose densities sum up into 
an azimuthally uniform distribution of the total intensity. 
However, multipole and necklace-ring vector solitons are 
subject to an azimuthal instability, except of the dipole- 
mode vector soliton [l9| and vortex-mode soliton close to 
a bifurcation |2^|. A particular class of such (generally, 
unstable) solutions is a vector vortex soliton with equal 
amplitude distributions in both components 
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A challenging issue is a possibility of the existence of 
stable vectorial solitons of the (to, —to) type, that would 
feature the same annular shape (with a hole in the cen- 
ter) as the bright scalar or vectorial vortices with the 
spins, respectively, m or (to, to), but with zero net spin. 
The possibility of the existence of such objects is ob- 
vious if the coupling between the components is of the 
XPM (cross-phase-modulation) type, i.e., insensitive to 
their relative phase - then, there is no difference in the 
shape between stationary vortex solitons of the (to, to) 
and (to, —to) types. However, in the simplest collapse- 
free model with saturable nonlincarity, which is based 
on XPM-coupled nonlinear Schrodingcr (NLS) equations 
for local slowly varying amplitudes £a,2, with terms like 
£1,2/ (l + |£i| 2 + I-E2I 2 ), the (m, — m) vortices are un- 
stable [as well as their (to, to) counterparts], although 
it was demonstrated that the instability may be partly 
suppressed [l4]. A stabilizing effect of the incoherent 
interaction of counter-rotating vortices was also demon- 
strated in an anisotropic photorefractive self-defocusing 
medium [l5|. 

Compound vortices with the spin components 
(1,-1,0) and (1,1,2) were also studied in the three- 
wave model of the type-II second-harmonic generation, 
with two components of the fundamental wave, and one 
second-harmonic component |l6tll7j . It was shown that, 
in this model per se, vortices of both types are unstable - 
against fusion into an ordinary zero-spin soliton, or split- 
ting, respectively. The addition of a stabilizing repulsive 
cubic interaction makes the life expectancy of the vor- 
tices much longer, but no case of complete stabilization 
has been found fl7| . 

In this work, we demonstrate that, in contrast with all 
the previously studied models, in the two-component CQ 
model the (m, — m) vortex solitons are rigorously stable 
in a certain parameter region. In fact, this result opens 
up a new class of stable 2D solitons with hidden vorticity 
(HV). 

We start with a general two-component CQ model 
in the rescaled form, which describes spatial evolution 
(along the propagation coordinate z) of the light beams 
in a bulk medium poj]. 

id z E h2 + AE 1:2 + (\E h2 \ 2 + a\E 2il \ 2 ) E h2 

- (\E h2 \ 4 + 2P\E 1 \ 2 \E 2 \ 2 + (3\E 2A \ A ) E l<2 = 0, (1) 

where E\^ 2 {x, y, z) are the local slowly varying ampli- 
tudes of the two waves, and the Laplacian A is the diffrac- 
tion operator acting on the transverse coordinates (x, y). 
The real parameters a and (3 are the cubic and quintic 
XPM coefficients, respectively. 

The CQ nonlinearity was experimentally observed in 
isotropic media, such as glasses |2l] and organic materi- 
als [23; these materials also feature nonlinear loss, but 
the consideration of the conservative model is justified, 
as the characteristic soliton's length can be made essen- 
tially smaller than the absorption length, or the loss may 
be compensated by gain. In that case, E\ and E 2 may be 



realized as orthogonally polarized waves, with a = 2/3 
for linear and a = 2 for circular polarizations, respec- 
tively. The latter case pertains as well to a set of two 
waves with different carrier wavelengths. Besides that, 
the model may serve as a crude isotropic approxima- 
tion for the description of photorefractive media in the 
low-saturation regime, which is characterized by equal 
strengths of the XPM and SPM (self-phase-modulation) 
nonlinearities, a = (3 = 1 2]. Thus, different values of a 
and [3 are physically relevant. 

It is well known that fundamental (m — 0) solitons, 
both scalar and vectorial [the latter ones, of the type 
(0, 0)], are stable in models with any saturable nonlinear- 
ity, including the CQ model (for the isotropic scalar sat- 
urable model, this was demonstrated more than a decade 
ago (23). Recently, stability of a vectorial soliton which 
is fundamental (m = 0) in one component, and carries 
the vorticity to = 1 in the other, i.e., a "vortex-mode 
soliton" of the (0, 1) type, has been pred icted close to a 
bifurcation in the saturable medium [24|. More complex 
vectorial solitons, combining a zero- vorticity configura- 
tion in one component and a topologically charged one 
in the other, are also possible in saturable media [Tl| . 
Such solutions definitely exist in the CQ model too, and 
they are plausibly stable in a broad parametric area. 

II. STATIONARY SOLUTIONS AND 
STABILITY ANALYSIS 

In this work we focus on the vectorial vortex solitons 
of the (m, ±m) types with symmetric components, which 
carry equal wave numbers k (non-symmetric solutions 
with ki k 2 are possible too, but we do not expect their 
properties to be drastically different from those of the 
symmetric solitons): 

(S)-"M-p(««)(5lSt)). < 2 > 

where r and <p are the polar coordinates in the plane 
(x,y), and the real function V obeys the equation 

kV = D m V+{l + a)V 3 -(l + 3f3)V 5 , (3) 

where D m = d 2 /dr 2 + r~ 1 d/dr — m 2 r~ 2 . Using the 
transformation V = R(l + a)^ 2 (l + 3/3)' 1 / 2 , r 
r(l + a) -1 (l + 3/3) 1 / 2 , and k — > fc(l + a)~ 2 (l + 3/3), 
we cast Eq. (0 in the form 

kR = D m R + R 3 — R 5 , (4) 

which is supplemented by the boundary condition R ~ 
r' m ' at r — ► 0. For r — > 00, there are two types 
of solutions to Eq. f4"|. coexisting in the medium with 
competing nonlinearities |25j| . viz., bright solitons with 

R ~ exp (—^/krj jsfr, and dark solitons with R 2 (r = 
00) = (1 + vl — 4fc)/2. Integral characteristics of the 
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FIG. 1: The wavenumber k (a) and Hamiltonian (b), calcu- 
lated as per Eq. (Hj, versus total power, P = Pi + P% % for the 
vectorial vortex solitons of the (m, ±m) types. 



bright vectorial soliton are its partial powers in both com- 
ponents, 



Pl,2 = 



2tt 



1 



rdrR 2 (r). 



(5) 



Global characteristic of the soliton families, in the form 
of dependencies k(P), are displayed in Fig. ^a) for 
a = f3 = 1, here the total power P = Pi + P 2 . The cutoff 
(largest possible) value of k for the bright-soliton fam- 
ily is the same as for the family of commonly known ID 
solitons in the CQ model, which is k (P = oo) = 3/16 = 
0.1875; at this value of k, the bright solitons become in- 
finitely broad (approaching a finite maximum amplitude, 
]? max = v3/2), i.e., they go over into dark solitons. We 
note that the Vakhitov-Kolokolov criterion, dk/dP > 0, 
which is a necessary condition for stability of the solitons 
[18|. is satisfied for all these solutions. Actually, it guar- 
antees the stability against radial perturbations, but not 
against symmetry-breaking ones, which are known to be 
most dangerous for the stability of vortex solitons . 

Note that the partial angular momenta of the compo- 
nents are not conserved independently, the corresponding 
integral of motion being the total angular momentum, 
AI = Mi(z) + M2(z), where the initial partial momenta 
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(recall tp is the angular coordinate). Two distinct values 
of total angular momentum, M = 2mP and M = 0, cor- 
respond to two types of solutions, (to, to) and (m, —to). 
Obviously, both solitons coincide in shape, therefore the 
diagrams for them shown in Fig.^are identical. However, 
stability of the two types of the solutions may be com- 
pletely different. An attempt to describe such a differ- 
ence was made, based on the "thin-ring" approximation, 
in a saturable medium [T3J; however, all the vortices are 
unstable in that medium. The "thin-ring" approxima- 
tion definitely does not apply to the "wide-ring" vortex 
solitons studied here. 

The most important results are adequately represented 
by the case of a — [3 — 1, on which we focus below. In 
this case, the Hamiltonian of Eqs. QJ is 



H = J (|V£i| 2 + |V£ 2 | 2 -±/ 2 + ±/ 3 )dr 



(7) 



with the total intensity I = \Ei\ 2 + |i?2| 2 - The transfor- 
mation of variables which leads to the normalized equa- 
tion J3J amounts to V(r) = R(r)/y/2, so that I = R 2 . 
With the latter substitution, Eq. Q reduces to the 
Hamiltonian of a scalar vortex with charge m, 



H = 2tt 



dR 

dr 



m -,2 1 n2 



-R z 
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2*" + §*" 



rdr. 



Perturbed vortex-soliton solutions are sought for in the 
form [cf. Eq. for the unperturbed ones] 



Ei = exp(ifcz + imtp)[R(r)/V2 + f + g*[ 
E 2 = exjp(ikz ± irrup) [R{r)/\2 + p + a* 



(9) 



here {f,g,p,q} = {f(r),g{r),p(r),q(r)}exp(X s z + isip) 
with complex eigenvalue X s and an arbitrary integer az- 
imuthal index s. Substitution of these expressions into 
linearized equations Q yields a system 
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(10) 



where L ± = D m±s - k + R 2 (3/2 - 2i? 2 ) and A = 
R 2 (1/2 - R 2 ). The signs ± in Eq. JTUJ) correspond to 
the two states (to, ±m). 

Note that in the case of (to, to) solutions, the matrix in 
Eqs. (|10fl has a block [2 x 2] structure, hence the eigen- 
modes degenerate (/ = p and g = q), and the linear 
stability problem reduces to the one for the scalar vor- 
tex, cf. Refs. 0,0- Therefore, the stability properties 
of the (to, to) vectorial vortices are completely identical 
to those of their scalar counterparts. However, the de- 
generacy does not take place for the HV solitons of the 
(to, —to) type, which clearly shows the difference in the 
stability problem for the two types of vectorial vortex 
solitons. 

Stability eigenvalues were found from numerical solu- 
tion of Eqs. i|10|) . In Fig. [21 we display dependencies 
of the eigenvalues with different values of the azimuthal 
index s on the wavenumber k for the vectorial vortex 
solitons of the (1,1) and (1,-1) types. The maximum 
growth rate is found for the modes with, respectively, 
s = 2 and s = 3. 

The above results comply with direct simulations of 
the evolution of the vortices shown in Fig. [3] for k = 0.1, 
when the linear stability analysis predicts that the soli- 
tons of both types are unstable. The symmetry-breaking 
instability modes of the (1,1) soliton in both components 
are identical, therefore in the Fig. [31(a) we display only 
one of them. The observed dynamics of two splinters 
(which are zero-vorticity vector solitons), generated from 
this soliton, is exactly the same as was observed for the 
scalar vortex: the splinters fly away in tangential direc- 
tions 
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FIG. 2: Growth rates of perturbation eigenmodes with dif- 
ferent values of the azimuthal index s (indicated next to the 
curves) for the vectorial vortex solitons of the types (1, ±1). 




FIG. 3: Instability-induced evolution of the vectorial vortex 
solitons with k — 0.1. (a) One of two identical components 
of the (+1, +1) soliton. (b) Two components of the (+1, —1) 
soliton with implicit vorticity. 
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FIG. 4: The same as in Fig.|5|for the vectorial vortex solitons 
of the types (m, ±m) for m = 2, 3. 



On the contrary, for the HV solitons [the (to, — m) type], 
the stability interval never reaches the cutoff value, i.e., 
dark vortices of the same type are always unstable. In 
fact, stability intervals for vortex solitons terminating at 
k smaller than the cutoff value ( "stability islands" ) have 
never been reported before. In the case of the higher- 
order solitons (with m = 2 and 3), the (to, —to) HV soli- 
tons have a smaller stability domain than their (m, to) 
counterparts. It is noteworthy too that the azimuthal 
index s of the most unstable eigenmode depends on the 
type of the soliton, and for the higher-order ones, of the 
types (2, ±2) and (3, ±3), the most dangerous value of s 
depends on k too. 



III. CHARGE FLIPPING 



A totally different scenario is observed in Fig. [3] (b) , 
where, in each component, three splinters of the initial 
HV soliton start to move in radial directions (cf. similar 
observations in Refs. [Hill). At this stage of the HV- 
soliton's break-up, as it seen in the panel corresponding 
to z — 180 in Fig. (b), the triangular sets of the splin- 
ters in the two components are slightly misaligned. With 
the further propagation, the separation of the splinters 
ceases, and they eventually fuse into a spinless (0, 0) vec- 
torial soliton. 

Figure 0] displays the instability growth rates found 
from Eq. (|10|) for higher-order vortex solitons, of the 
(to, ±to) types, for to = 2 and 3. As is seen from these 
figures, each type of the soliton has its stability area, 
as summarized in Table I. From these results, we con- 
clude that, for all the solitons with explicit vorticity [the 
(to, to) type], the stability regions extend up to the cut- 
off value, k — 0.1875, which implies that these solitons 
continuously carry over into stable vortices of the dark- 
soliton type, similar to what is known is the scalar case. 



Performing numerical analysis of the linear-stability 
problem, based on Eq. I|10|) . with higher accuracy we 
have found additional very small unstable eigenvalues 
for even azimuthal indices (s = 2to), with the growth 
rate ~ 5 x 10~ 4 , which were not visible in Figs. |2Jand0| 
(but which definitely exceed a possible numerical error). 
In Fig. |SJ we plot bifurcations which give rise to these 
miniscule eigenvalues (which are ordinary bifurcations of 
the pitchfork type), for the solitons with to = 1,2, and 
3. As a result of these additional very weak perturba- 
tions (which may be unobservable in the experiment), 
the corresponding rigorously defined stability regions for 
the (to, —to) HV solitons are considerably reduced with 
respect to their (to, to) counterparts, see Table I for a 
summary of the output of the stability calculations. 

To visualize the development of the weak instability 
revealed by Fig. \5\ we simulated the propagation of the 
solitons of both types (1,±1) with k = 0.14, adding 10% 
of the random noise. For the solitons with the explicit 
vorticity, the instability growth rate is ReA = 0.0373, 
and the dynamics follows the "usual" break-up scenario 
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FIG. 5: The growth rate of the small "internal" instabil- 
ity eigenmodes for the hidden-vorticity solitons of the type 
(m, — m) with m = 1,2, and 3, the corresponding azimuthal 
indices being s = 2m. 



shown in Fig.0](a), therefore we do not display it again. 
The growth rate of the same instability mode (with 
s — 2) but for the HV soliton of the (1, —1) type is two 
orders of magnitude smaller, ReA ~ 0.00036, therefore, 
noticeable development of the instability should be ex- 
pected after having passed the distance ~ 10 4 . Although 
so large propagation distances can hardly be achieved ex- 
perimentally (in the experiment, these solitons will seem 
as completely stable ones) , the issue is of principal inter- 
est, therefore we have completed the numerical analysis 
and found the soliton's dynamics of a novel kind. The 
results are summarized in Fig. where we show the in- 
tensity and phase distributions for both components up 
to z = 175000. 

To explain the complex internal dynamics observed in 
Fig. El we recall the theory developed for necklace-ring 
vector solitons in Ref. |l2(. The HV solutions that we 
consider here belong to a particular branch of a much 
broader class of necklace-ring vectorial solitons, namely, 
the ones with equal powers in both components, P\ — 
P-2 = P/2. In the most general case, these solutions can 
be represented as [cf. Eq. ©] 



Ei \ n 



' = i?(r)e lfc2 0*$ , . , 
E 2 J ysm (rrvp) 

Here the matrix 

= 



% Vl Q 



(11) 



(12) 



with arbitrary constants 0\ t2 indicates the phase invari- 
ance of the solutions, which is the symmetry property 
amenable for the conservation of the partial powers -Pi,2- 
The matrix <f> is simply a rotational transformation in 
the transverse plane, 



* = 



cos (mipo) sin (rrupo) 
-sin(m<^>o) cos(m<po) 



(13) 



where the arbitrary constant angle ipo reflects the rota- 
tional invariance responsible for the conservation of the 
total angular momentum. Note that, for radially sym- 
metric solutions Eq. @, the transformation ip — > ip— ipo, 
described by Eq. I|13fl . is equivalent to the phase shift 
accounted for by Eq. I|12|) with 9%^ = Tra^o. 
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FIG. 6: Evolution of the "internally unstable" vortex soliton 
with the hidden vorticity, of the (1, —1) type, for k = 0.14. 
The propagation distance is z = (a), 60000 (b), 120000 (c), 
and 175000 (d). The dynamical exchange of the spins 5*1,2 
(solid lines in (e)) is accompanied by charge flipping (sign 
reversal of the vorticity in each component), as clearly visible 
in the phase diagrams for arg(£i i 2) in (a)-(d). The total 
spin S = (Si + S2)/2 keeps its zero value (the dashed line 
in (e)). The net intensity, |i?i| 2 + I-E2I 2 , keeps its azimuthal 
homogeneity in the course of the entire propagation, as shown 
in (f) for z — 175000. Vertical dotted lines correspond to 
snapshots (b) and (c). 



With the accuracy of arbitrary phase shift being al- 
ready absorbed by the matrix 0, the linear transforma- 
tion \& in Eq. (|11|) is given by 



cos ip i sin ip 
sin ip ±z cos ip 



(14) 



where the constant parameter ip describes the rotation 
in the space of the components {E\,E 2 ), similar to the 
Manakov system, and corresponds to the conservation if 
the "isotopic spin", i J {E\E 2 — ElE 2 }dr. For any value 
of ip, the latter expression is zero in our case. In addi- 
tion, the parameter ip uniquely defines the initial values 
of the partial spins Si = Mi j Pi [recall the partial angular 
momenta were denned in Eq. JfjJ]: 



Si = msin(2-0), S 2 = ±Si, 



(15) 



and, therefore, it determines the total spin (dynamical 
invariant), S = M/P = |(Si + S 2 ) = |(m ± m) sin(2-0). 
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Among possible stationary solutions conforming to 
Eq. (|llfl are those with zero, fractional (0 < S < to), 
and integer total spin (S — to). We focus here on two 
cases which correspond to Eq. @: the (to, to) type of the 
solutions, with the total angular momentum attaining its 
maximum possible value, M = mP (i.e., S = to), and the 
HV solutions of (to, —to) type, with M = S = 0. Both of 
them represent radially symmetric vector vortices with 
\ip\ = tt/4; the transformations ip — > —ip, to — > —to, and 
(to, ±to) — > (— to,=f?ti) are all equivalent. 

Solutions with the explicit vorticity, S = to, corre- 
spond to the upper sign in Eq. ljbf|) • The partial spins as- 
sume the maximum possible values in this case, Si^ = to, 
thus the only corresponding configuration is the one with 
\ip\ = 7r/4, and the vortex soliton of this type always has 
axially symmetric (ring-shaped) components. In addi- 
tion, the exchange of angular momentum between com- 
ponents is forbidden in this case. 

The HV solutions with 5 = are drastically differ- 
ent. They correspond to the lower sign in Eq. I|14|) for 
an arbitrary value of ip. These include the HV vectorial 
solitons for ip = ±7r/4, and also solutions with inten- 
sity distributions in the two components in the form of 
two crossed multipoles, for ip = and = 0. Similar 
solution with to = 1, or a dipole-dipole vectorial soli- 
ton, was investigated theoretically and experimentally 
in Refs. [T3 . Il3| and found to be azimuthally unstable 
in saturable medium. The ^-values from the interval 
— 7r/4 < ip < 7r/4 determine the depth of the azimuthal 
modulation in each of the two components, that sum up 
into the azimuthally uniform distribution of the total in- 
tensity I = R 2 12] . Thus, the continuous soliton family 
includes crossed multipoles with different values of the az- 
imuthal modulation depth and opposite fractional values 
of the partial spins. It is important to note that, because 
all the dynamical invariants do not depend on ip, for the 
whole branch of the HV vectorial solitons, parametrized 
by ip, the dynamical exchange of the angular momentum 
between components is possible. 

Using the decomposition in Eq. Qllfl. we can now ex- 
plain the dynamics in Fig. 03 in terms of a slow evolution 
of the parameters tpo and ip, as shown in Fig. [7] In other 
words, instead of the modulational instability along the 
azimuthal direction, that would lead to fragmentation 
of the initial annular soliton, the instability modes from 
Fig. El initiate slow sliding of the solution across the con- 
tinuous manifold with the independent parameters ipo 
and ip. At each stage of the evolution, such as those cor- 
responding to the frames (a)-(d) in Fig. 03 we observe 
a slightly perturbed stationary solution with varying ipo 
and ip (the arbitrary phases are of no importance 
because the XPM interaction between the components is 
phase-insensitive) . 

Indeed, the random noise on the level of 10% of the 
soliton's amplitude, added to the HV soliton in Fig.EJa), 
quickly dissipates, and, in the course of the first several 
thousands units, the vector HV soliton propagates with- 
out any noticeable change. Then, as is seen in Fig. Ufa), 
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FIG. 7: Evolution of the parameters ip from Eq. 1141 in (a) 
and ipo from Eq. I13H in (b), corresponding to the adiabatic 
"sliding" across the soliton family, denned by Eq. (li lt in the 
course of the evolution displayed in Fig.|5] In (c), the angular 
velocity of the internal rotation is shown; strong numerical 
noise on this curve reflects the uncertainty of the value of the 
parameter ipo (which is arbitrary, in the general case). Ver- 
tical dashed lines indicate the charge-flipping (spin-reversal) 
points, close to those shown in Figs. 03(b) and (d). 



the parameter ip decreases and the components assume 
a shape of two crossed dipoles. When the modulation 
depth reaches its maximum for ip = 0, the solution is, 
simply, {E 1 ,E 2 } = R(r){cos((p-tp Q ),sin((p-(p )}, hence, 
at this point, neither component contains any vorticity, 
as is indeed seen in Fig. |SJb). Because the parameter tpa 
is arbitrary, its particular value at this stage (Fig. 0(b)) 
depends on the initial noise. 

With the further propagation, the components £1.2 
almost restore their initial annular shape and develop 
the phase dislocations corresponding to the vortex cores 
afresh, which are opposite to initial ones, cf. Figs. 03(a) 
and (c). This phenomenon, the "charge flipping", was 
recently predicted to occur in a significantly different sys - 
tem, namely, vortices in nonlinear photonic lattices j26j. 
The latter system does not conserve the angular momen- 
tum at all, because the rotational symmetry is broken 
by the lattice. Nevertheless, the similarity with that sys- 
tem, which seems to be important for the effect to occur, 
is the presence of two interacting sub-systems which are 
given the freedom to exchange the angular momentum: 
the vortex and the lattice in Ref. j2(|, and the two com- 
ponents of the vortex in the present case. It is commonly 
known from the studies of vortices in linear optics |27| 
that vortices in the phase fronts can annihilate or be born 
in pairs. The charge-flipping phenomenon introduces a 
new mechanism of such transformation in the nonlinear 
setting, through the exchange of the angular momentum 
between two nonlinear ly (XPM) coupled subsystems. 

After the first charge flip, the components do not fully 
restore their annular shape (Fig. 03 (c)), and the maxi- 
mum value of the partial spins which have the opposite 
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FIG. 8: Final stage with soliton "explosion" after quasi-stable 
evolution in Fig. 0for the propagation distances z = 179500 
in (a), 179600 in (b), and 179700 in (c). 



signs, S% = —Si w 0.83 with ip w —0.462, is attained at 
z ps 117300 (Fig. 0^a)). At the same time, the dipoles 
in both components start to rotate slowly (the angle ipo 
increases, see Fig. 0(b)). We can introduce, therefore, 
the angular velocity of the rotation, depo/dz, as shown in 
Fig. 0c). However small, it demonstrates an important 
feature of the correlation between the internal "degrees 
of freedom" ip and ipo . Indeed, for the exact stationary 
solutions, these two parameters are independent, while 
for the perturbed solutions in Figs.0and0 they become 
coupled through the growing instability modes. In par- 
ticular, points where the angular velocity vanishes corre- 
spond to dipole-dipole soliton, with zero vorticity in both 
components, while fastest rotation is achieved when the 
partial spins in the components attain maximum abso- 
lute values. 

We continued the simulations, and eventually observed 
a break-up of the vortex soliton, as shown in Fig. EI It oc- 
curs within the distance of several hundreds of the propa- 
gation units, i.e. three orders of magnitude smaller than 
the previous stable propagation, thus it may be consid- 
ered as an "explosion" . We have checked the evolution of 
the integral of motions at this stage, and, in particular, 
observed perfect conservation of the total spin, Fig.0(d), 
which rules out a numerical error as a probable cause of 
the "explosion". 



(a) 



(b) 



Power, P 



Hamiltonian, H 




FIG. 9: Evolution of the total integral power (a), Hamiltonian 
(b), and accordingly defined propagation constant (c) in the 
course of the propagation shown in Figs. and 



In Fig. we plot the total integral power and Hamil- 
tonian, calculated in the course of the propagation, and 
notice small changes which naturally occurs due to the 
radiation emission from perturbed soliton. The breakup 
is accompanied by strong radiation which leads to sharp 
changes in the final segment of the diagrams. Using the 
relations displayed in Fig. we restore the correspond- 
ing value of the propagation constant k and plot it in 
the Fig. (c). As one can see, the propagation constant 
decreases, parallel to the power loss due to the radia- 
tion. The explosion occurs when the propagation con- 
stant reaches the value k ps 0.13627, and the splitting 
follows the scenario already observed in Fig. - three 
splinters fly away along radial directions. We conclude 
that the sudden splitting of the otherwise "intrinsically" 
unstable HV soliton happens because, as a result of the 
slow evolution of its integral characteristics due to the 
continuing radiation loss, it hits a boundary of the "ex- 
ternal" instability domain, after which it breaks apart. 
The latter boundary is found from the linear-stability 
analysis to be at k = 0.13582 (see Table I), and we stress 
remarkable agreement and accuracy of the numerical pro- 
cedure: both values coincide up to 10~ 3 , k « 0.136. 

Deeper in the linear-stability domain, for instance for 
k = 0.16, the HV soliton of the type (1, —1) demonstrates 
no sign of instability for any reasonable propagation dis- 
tance despite addition of initial noise, because the cor- 
responding instability mode with s = 2 in Fig. has 
a vanishingly small growth rate. We reflect the fact of 
small "internal" instability by excluding corresponding 
domain from the final summary in Table I. 

We believe that the internal dynamics of the vectorial 
soliton reported above (Figs. and for the particular 
case of m = 1 and the CQ nonlinearity can manifest 
itself as a generic effect, for higher topological charges 
and in other systems, for example in the mixture of Bose- 
Einstein condensates 1281. 
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FIG. 10: (a) Real (solid) and imaginary (dashed) parts of the 
instability eigenvalues with the maximum growth rate for the 
vectorial solitons of the (1, ±1) types, (b) An additional bi- 
furcation (magnified in the inset A) of the perturbation mode 
with s = 1, which occurs for the hidden- vorticity soliton of 
the (+1,-1) type close to the cutoff. This instability gives 
rise to a shift of the soliton's core off the center. 



IV. INSTABILITY OF THE VORTEX CORE 
CLOSE TO CUT-OFF 

As it was already stressed, Figs. [5] and ^demonstrate 
that, in contrast to the stability domain for the solitons of 
the (m, m) type, the stability region for the HV solitons 
of the (to, — m ) type does not extend to the cutoff point. 
In this section, we aim to study the instability of the 
HV solitons close to this point. The simulations demon- 
strate that the corresponding instability mode with the 
azimuthal index s = 1 leads to a shift of the vortex core. 
A possibility of this specific instability was earlier studied 
in Ref. |29j for vortices in scalar models, using an ana- 
lytical approximation for very broad annular solitons. It 
was concluded that such an instability may occur, but no 
linear unstable mode corresponding to the core shift was 
found [as before, the scalar case is exactly tantamount to 
the symmetric vectorial vortices of (to, to) type consid- 
ered here]. In contrast to that, Fig. I1UI shows an explicit 
example of such an instability mode for the HV solitons 
of the (1, —1) type. 

Numerical development of this instability is displayed 
for k = 0.18 in Fig. ^2 Adding 20% of an initial random 
noise (Fig. El ( a )) does not strongly affect the dynam- 
ics — for up to 4000 propagation units, it shows no sign 
of instability. The only visible action of the perturba- 
tion is excitation of the internal modes of the vectorial 
soliton, corresponding to the purely imaginary eigenval- 
ues in the linear-perturbation spectrum (such eigenval- 
ues are shown in Fig. 1101 (a)). The pattern displayed in 
Fig. 1111 (b) periodically repeats itself during the prop- 
agation. Similar long-lived internal modes have been 
recently observed in perturbed evolution of scalar CQ 
vortex solitons jsjj. It is noteworthy too that, for the 
symmetric vectorial vortex shown in Fig. 1111 the total 
intensity shows no sign of azimuthal modulation, as the 
intensity of the components sum up to the axially uni- 
form distribution, similar to the case of "intrinsic" insta- 
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FIG. 11: Dynamics of the hidden-vorticity soliton of the 
(1,-1) type with k = 0.18 and initial 20% random-noise per- 
turbation. The propagation distance is z — in (a), 500 in 
(b), 4300 in (c), and 5000 in (d). As before, the frame (e) 
demonstrates the evolution of the partial spins (solid lines) 
in the two components and conservation of the total angular 
momentum (dashed line). 



bility shown in Fig. |SJ It is expected that the full set 
of internal modes of the vectorial vortex solitons should 
inherit all symmetry properties of the underlying station- 
ary solutions; however, no verification of this assumption 
has been reported thus far. 

At the distance z — 4300, the growth of the unsta- 
ble mode results in the separation of the dislocations in 
the two components, see Fig.^2( c )- The total- intensity 
distribution remains unmodulated, and the vortex core 
is invisible, because it is covered, in the total-intensity 
distribution, by the mutual displacement of the compo- 
nents: actually, an intensity maximum in one component 
lies on top of a minimum in the other. The shift of the 
vortex core quickly increases, and it moves to the outer 
edge of the soliton. As a whole, the vectorial soliton re- 
mains localized, as is seen in Fig.llll(d). and it possesses 
no vorticity, as the spin diagrams in Fig. 1111 (e) demon- 
strate. Strong modulation of the components inside the 
soliton persists for the a long propagation distance after 
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TABLE I: Stability intervals and their size relative to the 
overall existence interval, < k < 3/16 = 0.1875, for the 
vectorial vortex solitons of sundry types. 





unstable 


stable 


%% 


(2,2) 
(3,3) 


< k < 0.14855 
< k < 0.16190 
< k < 0.17005 


0.14855 < k < 0.18750 
0.16190 < k < 0.18750 
0.17005 < k < 0.18750 


20.8% 
13.7% 
9.3% 


(1,-1) 
(2,-2) 
(3,-3) 


< k < 0.13582 
< k < 0.14884 
< k < 0.15866 


0.16163 < k < 0.17945 
0.15620 < k < 0.15940 
0.15866 < k < 0.15973 


9.5% 
1.7% 
0.57% 



the vortex annihilation. We conclude that, in this case, 
the result of the development of the instability mode is 
establishment of a breather, i.e., a strongly perturbed 
and oscillating zero-spin soliton. 

V. CONCLUSION 

We have demonstrated, for the first time, that 2D spa- 
tial solitons of the annular shape, carrying zero total vor- 
ticity, may be stable in the CQ medium, being supported 
by the hidden (implicit) vorticity, in quite a broad region. 
The output of the linear-stability analysis is summarized 
in Table I, where the stability domain and its size rela- 
tive to the existence domain are shown. There is a single 
border between instability and stability regions for the 
cxplicit-vorticity solitons, of the (to, to) type. The stabil- 
ity domain for this type of the vectorial solitons extends 
up to the point of the transition to dark-soliton vortices 
(cutoff). For the implicit-vorticity solutions of (to, — m) 
type, the situation is more complex. We identify the re- 



gion of the relatively strong "external" instability (shown 
in Table I), where unstable vectorial solitons with both 
implicit and explicit vorticity split into a set of fragments, 
the number of which is equal to the azimuthal index of 
the fastest growing mode of small perturbations. For 
larger k, and only for the hidden- vorticity solitons of the 
(to, —to) type, there exist a domain of the very weak "in- 
trinsic" instability, where the vectorial soliton as a whole 
remains robust, while its components undergo a very slow 
internal evolution, "sliding" through the solution family 
and exhibiting exchange of the angular momentum and 
charge flipping. In a vicinity of the cutoff, where the 
explicit vortices are stable, the hidden-vorticity solitons 
reveal a weak-instability mode that results in splitting of 
the phase dislocations in the two components, somewhat 
similar to the splitting of multiple-charged dark vortex 
solitons in the scalar model. Therefore, the stable soli- 
tons with hidden vorticity may be regarded as "excep- 
tionally bright" objects, unlike the familiar solitons car- 
rying explicit vorticity, which always have stable dark- 
vortex counterparts. With the increase of the integer 
vorticity to, the stability regions of the vectorial solitons 
with both explicit and implicit vorticity quickly shrink. 
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